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ABSTRACT 

Context. A promising source of the positrons that contribute through annihilation to the diffuse Galactic 511 keV emission is the 
y8*-decay of unstable nuclei like '*Ni and '''^Ti synthesised by massive stars and supemovae. Although a large fraction of these 
positrons annihilate in the ejecta of SNe/SNRs, no point-source of annihilation radiation appears in the INTEGRAL/SPI map of the 
511 keV emission. 

Aims. We exploit the absence of detectable annihilation emission from young local SNe/SNRs to derive constraints on the transport 
of MeV positrons inside SN/SNR ejecta and their escape into the CSM/ISM, both aspects being crucial to the understanding of the 
observed Galactic 511 keV emission. 

/Wefhods. We simulated 511 keV lightcurves resulting from the annihilation of the decay positrons of ^^^Ni and ■*''Ti in SNe/SNRs and 
their surroundings using a simple model. We computed specific 511 keV lightcurves for Cas A, Tycho, Kepler, SN1006, Gl.9+0.3 
and SN1987A, and compared these to the upper-limits derived from INTEGRAL/SPI observations. 

Results. The predicted 511 keV signals from positrons annihilating in the ejecta are below the sensitivity of the SPI instrument by 
several orders of magnitude, but the predicted 511 keV signals for positrons escaping the ejecta and annihilating in the surrounding 
medium allowed to derive upper-limits on the positron escape fraction of ~13% for Cas A, ~12% for Tycho, ~30% for Kepler and 
-33% for SN1006. 

Conclusions. The transport of ~MeV positrons inside SNe/SNRs cannot be constrained from current observations of the 511 keV 
emission from these objects, but the limits obtained on their escape fraction are consistent with a nucleosynthesis origin of the 
positrons that give rise to the diffuse Galactic 511 keV emission. 

Key words. Gamma rays: general - ISM: supernova remnants - Nuclear reactions, nucleosynthesis, abundances - Astroparticle 
physics 



1. Introduction 



One of the most valuable contributions of the INTEGRAL 
space-borne gamma-ray observatory so far is the cartography 
of the 511 keV galactic emission from electron-positron anni- 
hilation (^SdTsedeFeFarjSOOS Weidenspointner et al.|2008a i. 
This achievement was obtained from several years of observa- 
tion with the high-resolution SPI spectrometer and provides us 
with the most accurate view on leptonic antimatter in the Galaxy. 
This result, and its expected improvement over the coming years 
as exposure grows (INTEGRAL is planned to operate till at least 
2012), will very likely help clarifying the origin of the positrons 
that annihilate in our Galaxy, an issue that has been lasting for 
several decades and that has proven over the years to be of inter- 
est for a number of different science fields, from astrophysics to 
particle and plasma physics. 

A review of the history and stakes of annihilation gamma-ray 
observations, from early balloon experiments to the SMM/GRS 
and CGRO/OSSE satellites, can be found in [Prantzos et al.| 
( 2010| l and Diehl & Leising (2009). The above-mentioned 
INTEGRAL result confirms the basic description of the emission 
that followed from the CGRO mission, a bright bulge typically 
10° in width on top of a faint disk. The positive latitude enhance- 
ment found from OSSE data is however not observed by SPI 
and instead, a longitudinal asymmetry in the inner disk is found 



( Weidenspointner et al.|2008a|b i. The spectrometric capabilities 
of the SPI instrument also allowed to derive a high-resolution an- 
nihilation spectrum, and the observed line profile indicates that 
most of the annihilation occurs through positronium formation 
and takes place in the warm phases of the interstellar medium 
Pean et al.|2006l ). 

The observed intensity distribution, with a high inferred 
bulge-to-disk ratio of 2-6 in luminosity (Weidenspointner et al. 
|2008bl l, is intriguing in that it does not correlate with any distri- 
bution of classical sources or potential annihilation medium. For 
instance, thermonuclear supemovae (SNIa), that are promising 
positron sources, are ~ 10 times more numerous in the disk than 
in the bulge and would give a bulge-to-disk luminosity ratio of 
~ 0.1 if the positrons were to annihilate not too far from their 
sources (Prantzos 2006 ). The peculiar morphology of the galac- 
tic 5 1 1 ke V emission has been invoked as an evidence for an ex- 
otic origin of the positrons, such as dark matter decay in the in- 
ner regions of the Galaxy ( |Boehm et al.|2004| . Another explana- 
tion put forward is the large-scale transport of positrons created 
by more classical sources (P rantzos 2006; Higdon et al.|2009 i. 
The latter alternative is supported to some extent by the spectro- 
scopic properties of the annihilation line. Indeed, most plausible 
positron sources are found in the hot phase, which occupies most 
of the interstellar volume, while most positrons are observed to 
annihilate in the warm phase ( Jean et aL]|2006 Churazov et al. 
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|2005| l; positron transport therefore seems to be required to a cer- 
tain degree. Overall, the discussions raised by the INTEGRAL 
observations illustrate that understanding the observed annihila- 
tion emission is both a problem of sources and transport of the 
positrons. 

Among the most plausible sources of galactic positrons is 
the yS^ decay of unstable nuclei synthesised by massive stars and 
supernovae, mainly ^^Ni, '^'^Ti and ^^Al. The existence and de- 
cay of the above-mentioned isotopes are directly established by 
gamma-ray l ine observ ations (I Leising & Share 1990[ Renaud 



cascade of MHD fluctuations is quenched at too large spatial 
scales as a result of ion-neutral collisions ( Higdon et al.||20d9 



let al.|2006bi [Martin el al. 2009b)|, and die observed or inferred 
isotope yields are fairly sufficient to account for the steady state 
annihilation rate of a few 10"^^ e^ s ' deduced from INTEGRAL 
observations. Under the assumption that positrons annihilate 
close to their sources, the estimated annihilation rate in the disk 
of the Galaxy agrees with the estimated pos itron production 
rate from ^^Al decay (Knodlse der et al.||2005 i. Because of the 
high bulge-to-disk ratio in luminosity, however, the origin of the 
positrons annihilating in the bulge of the Galaxy is much more 
uncertain. From COMPTEL observations, the initial energy of 
these positrons was constrained to be below a few MeV, other- 
wise they would give rise to a difflise MeV emission of in-flight 
annihilatiorQ (see for instance Beacom & Yiiksel 2006 Sizun 



et al.|2006) l. In that respect, the positrons from nucleosynthesis 



are good candidates since they are released with initial kinetic 
energies of ~ 1 MeV (but see Chernyshov et al.|2010[ ). 

In order to contribute to the 511 keV galactic annihilation 
emission, the decay positrons must first be slowed down to be- 
low ~ 100 eV and then annihilate in optically thin environments. 
Since the interstellar medium is fully transparent to 5 1 1 ke V ra- 
diation, this means that the parent radio-isotopes and/or their 
positrons should escape the dense stellar interiors where they 
have been produced. The slowing-down and annihilation of de- 
cay positrons may take a certain time (see |2.1| i, during which 
the initially relativistic particles can diffuse. The ultimate fate 
of decay positrons therefore depends on the characteristics of 
their transport, first in the expanding stellar ejecta of super- 
novae (SNe) and supernova remnants (SNRs) and then in the cir- 
cumstellar medium (CSM) and interstellar medium (ISM). The 
transport of MeV positrons, however, is as yet largely unsolved. 
The nature of their diffusion regime, collisionless (positrons ef- 
ficiently scatter off MHD waves) or collisional (positrons follow 
unperturbed trajectories along magnetic field line^, is still un- 
known. This can affect the characteristic path lengths travelled 
by positrons before annihilation by several orders of magnitude, 
from several pc to a few ten kpc. The problem stems firstly from 
our incomplete understanding of how MHD turbulence scatters 
particles through their interaction with MHD waves. Secondly, 
we have a poor knowledge of the MHD turbulence at small spa- 
tial scales (of the order of the Larmor radius of MeV positrons 
in ~ juG fields), where the decay positrons could resonantly in- 
teract with perturbation modes. Recent works indicate that the 



' Very energetic positrons can also give rise to a GeV Galactic dif- 
fuse emission through inverse-Compton and bremsstrahlung radiation, 
as do the ~ 10 times less numerous secondary positrons coming from 
hadronic interactions of cosmic-rays with ambient gas (see |Strong et al.| 
I004}|Aharoman & Atoyan|2000[ >. 



^ The term "unperturbed" refers here to the absence of efficient scat- 
tering, but the particle constantly loses energy along its trajectory. 
Coulomb or ionisation/excitation interactions involve relatively small 
energy variations for MeV positrons and act like continuous energy- 
loss processes with small pitch-angle scattering. The helicoidal motion 
of the particle around a magnetic field line is thus only marginally af- 
fected. 



Jean et al.|2009 |l. In this case, no collisionless diffusion through 
resonant interactions takes place and positrons can travel over 
large linear distances. Yet, many aspects of the issue should still 
be investigated, like non-resonant interactions with fast magne- 
tosonic modes, which could lead to efficient diffusion even if no 
magnetic fluctuations exist at the resonant scale, or perpendic- 
ular diffusion by field line wandering ( Prantzos et al.|[20T0 l. A 
complete description of positron transport will also need to take 
into account diffusion in momentum space (reacceleration). On 
the whole, the issue is far from settled and theoretical as well as 
observational efforts are still needed. 

A nucleosynthesis origin for the Galactic positrons therefore 
requires a sufficient fraction of decay positrons to escape the pro- 
duction sites of their parent isotopes and feed the galactic popu- 
lation, and a global transport of these particles in the ISM so that 
the distribution of their annihilation sites can account for the ob- 
served morphology. 

As we will see later, a large fraction of the positrons created 
in the decay of ^^Ni and '^"^Ti are thought to annihilate inside 
SNe and SNRs. Despite this fact, no point-like emission appears 
in the INTEGRAL 511 keV map. This non-detection may there- 
fore be translated into some evidence concerning the intertwined 
aspects of MeV positron transport and escape, both topics being 
crucial to the understanding of the observed Galactic annihila- 
tion emission. 



2. Decay positrons in SNe and SNRs 

Positron production by the three above-mentioned radio- 
isotopes occurs in the following way: 



56Ni 56, 



44i 



Ti 



85 yr 



Co > 

6h 



44 



Sc 



^''Fe-He+ (18.1%) 
'*'*Ca + e+ (94.3%) 



where the decay timescales and yS-fractions are given for each 
decay chain. The electron neutrinos and antineutrinos (required 
for conservation of leptonic number) and the gamma photons 
(from the de-excitation of daughter nuclei) have been omitted. 
The positrons from ^''Ni and '*''Ti are actually released in the de- 
cay of their daughter nuclei, but in the following we will keep 
the term ^''Ni and '^'^Ti positrons to emphasise the connection 
with nucleosynthesis products. 

In the context of SNe/SNRs, the specificities of these most 
likely contributors to the Galactic positron budget are the fol- 
lowing: 

1 . 5^Ni is synthesised by explosive Si-burning deep in the stel- 
lar ejecta during core-collapse and thermonuclear SN ex- 
plosions, with typical yields of a few 10"' Mq. It is short- 
lived and escape is therefore a quite critical problem as the 
dense medium of the late SN/early SNR is likely to stop all 
positrons and absorb some of the annihilation radiation. 

2. '^''Ti is synthesised by explosive Si-buming deep in the stel- 
lar ejecta during core-collapse and thermonuclear SN ex- 
plosions, with typical yields of a few 10^^-10"^ Mq. It is 
medium-lived and depending on the conditions, the escape 
fraction can vary between two extremes (no or total escape). 
Yet, the diluted ejecta at typical decay time is mostly trans- 
parent to gamma-rays and so an annihilation signal can be 
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expected from those positrons that cannot escape the stellar 
envelope. 

3. ^*A1 is synthesised at various stages of the evolution of a 
massive star with typical yields of a few lO^'^-lO"^ M©. It is 
ejected in the ISM through stellar winds and SN explosions. 
Since it is long-lived, escape is not a problem and all decay 
positrons are directly injected in the ISM. 

In the present work, we are interested in point-sources of anni- 
hilation radiation so we will not consider ^''Al positrons, which 
are expected to give rise to a diffuse annihilation signal (in the 
Galactic disk in particular). We therefore focus on ^^Ni and '^"^Ti 
positrons only in the context of SNe and SNRs. Before address- 
ing the specific case of positron transport and annihilation in 
SNe/SNRs, however, we describe in more detail the main phys- 
ical processes that rule the life of a MeV positron until its anni- 
hilation. 

2.1. Life and death of MeV positrons 

Once created in the yS-decay of a given radionucleus, positrons 
are slowed down by ionisation/excitation and Coulomb losses. 
If their kinetic energy is brought down to ~ 100 eV, they can rip 
electrons off atoms (a process called charge exchange) and form 
positronium, which is a short-lived electron-positron bound state 
(the overall process is termed "positronium formation in-flight"). 
Depending on spin orientation, positronium then quickly decays 
in two 511 keV photons (25% of cases) or three continuum pho- 
tons (75% of cases, with the total radiated energy being equal 
to 2 nieC^ in the positronium rest frame). The positrons that do 
not form positronium in-flight eventually thermalise and anni- 
hilate directly with bound or free electrons, thereby giving two 
511 keV photons, or through positronium formation by radia- 
tive recombination with free electrons (for a review of the vari- 
ous annihilation channels, see Guessoum et al. 2005). It should 
be noted here that the cross-sections for positronium formation 
by charge exchange are several orders of magnitude larger than 
for any other process. So unless the medium is strongly ionised 
(and hence charge exchange cannot efficiently take place), most 
positrons will form positronium and then quickly annihilate once 
they are slowed down to a kinetic energy of ~ 100 eV. If the 
medium is strongly ionised, positronium formation by radiative 
recombination of thermalised particles is the dominant process, 
provided the temperature is moderate (for temperatures above ~ 
10^ K, the Maxwellian distribution of particles is shifted to an 
energy range where direct annihilation with free electrons takes 
over). 

The physical properties of the medium therefore influence 
the channel by which most positrons annihilate, and hence the 
fraction of positrons that annihilate through 511 keV emission 
as opposed to continuum emission. Different channels also im- 
ply different lifetimes for the positrons. The typical time to 
slow down a MeV positron down to ~ 100 eV through ionisa- 
tion/excitation or Coulomb losses is ~ 10"^"^ n"' yrs (where n 
is the medium density in cm"^) and the time for complete ther- 
malisation is similar since the energy-loss rate strongly increases 
as the particle energy decreases. Then, the time for positronium 
formation by charge exchange is comparatively negligible so for 
this annihilation channel, the lifetime of positrons is imposed by 
the slowing-down time. In contrast, the time for direct annihi- 
lation of thermalised positrons with free or bound electrons is 
~ 10''"^ n"' yrs, so in the latter scenario the lifetime of positrons 
is imposed by the annihilation time. In all cases, however, the 
range of positrons is imposed by the slowing-down period. 



2.2. MeV positrons in SNe/SNRs 

A SN ejecta is a relatively dense medium in which the decay 
positrons will be slowed down quite efficiently in early times and 
then on longer and longer timescales as the stellar material gets 
diluted. For the simple case of a uniform ejecta in homologous 
expansion and under the assumption that positrons remain con- 
fined to the stellar ejecta and do not escape, the slowing-down 
time for a MeV positron injected at a given time after the explo- 
sion can be approximateqjby: 

tsd^lQ-'x^xP yrs (1) 

where M^j and E^j are the ejecta mass and energy respectively 
(in Mo and 10^' erg) and f the time since the explosion at which 
the positron is injected (in yrs). 

Owing to its very rapid expansion, the stellar ejecta is 
thought to be quite cold and hence neutral, which means that 
positrons are slowed down by ionisation/excitation losses and 
then annihilate mostly through positronium formation in-flight 
just after the slowing-down period. From the above formula, we 
expect the slowing-down time to be below a year during the 
first years after the SN explosion. Over this period, the annihi- 
lation lightcurve therefore follows the radioactive decay curve. 
Then, between 10 and 100 yrs, the slowing-down time increases 
up to few centuries and so the positrons created in this inter- 
val and trapped in the ejecta will experience a delayed anni- 
hilation over a characteristic timescale that corresponds to the 
typical ages of the youngest Galactic SNRs (a few centuries to 
a millenary). After a century, the slowing-down times become 
too long for annihilation to occur inside the SNR and the cor- 
responding positrons very likely escape. For the simple case we 
are considering here, we therefore expect most positrons from 
the short-lived ^^Ni to annihilate soon after their release in the 
ejecta. Yet, a tiny fraction of them will give rise to a delayed an- 
nihilation signal which, given the high yields of ''^Ni, may be in 
reach of modem gamma-ray instruments. "^'^Ti holds more poten- 
tial for emission from SNRs because of its ~ 100 yr lifetime. At 
early times, positrons annihilate nearly instantly, as in the case 
of ^^Ni, but the sizeable fraction of positrons created after a few 
decades will lead to a delayed annihilation, potentially observ- 
able today in the youngest Galactic SNRs. 

This scenario is rather simplified and many aspects of the 
problem could alter the above conclusions. First of all, the 
above results are valid for 1 MeV positrons whereas most decay 
positrons actually have lower energies of a few 100 keV. Then, 
we made the assumption that decay positrons are "trapped" in 
the ejecta, which appears as the most favourable case for efficient 
slowing-down and subsequent annihilation, but positron escape 
very likely occurs (and actually should occur if the Galactic an- 
nihilation emission, especially from the Galactic bulge, is to be 
traced to a nucleosynthesis origin), which would lead to reduced 
511 keV fluxes. On the other hand, escaping positrons that lost 
some fraction of their kinetic energy during ejecta crossing can 
achieve their slowing-down and annihilate in the circumstellar 
environment of the SNR. 

Our order-of-magnitude evaluation, together with the above 

^ The formula gives the slowing-down time assuming that the 
positron injected at a given time experiences the same ejecta density 
during its slowing-down. As time goes by and slowing-down time in- 
creases, this assumption becomes less and less valid and the estimated 
value turns to a lower-limit. 
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considerations, show that there is potential for detectable an- 
nihilation radiation from SNe/SNRs up to a few centuries 
after the explosion, but such an emission is not observed. 
INTEGRAL/SPI observations provide the most stringent con- 
straints on annihilation in SNe and SNRs so far, and there- 
fore call for a deeper interpretation of the results. In this 
work, we focused on the 511 keV emission from a few young 
and well-characterized SNRs to constrain the transport of 
positrons in the ejecta and in the CSM/ISM. Upper-limits on the 
511 keV annihilation flux from six young SNRs, Cassiopeia A, 
Tycho, Kepler SN1006, SN1987A and G1.9-H0.3, were derived 
from INTEGRAL/SPI data for various annihilation scenarios. 
These constraints were then compared to simulated annihilation 
lightcurves computed in the frame of simple descriptions for the 
SN/SNR evolution and CSM/ISM environment, and under vari- 
ous prescriptions for the transport of the particles. 

3. Past works 

In order to evaluate the contribution of decay positrons to the 
galactic annihilation radiation observed by the GRIS and OSSE 
instruments, Chan & Lingenfelter] ( |1993 1 performed a theoreti- 
cal study of the escape of ^°Ni and '*'*Ti positrons from SNe of 
both types (thermonuclear SNIa and core-collapse SNe). Their 
results indicate that between 0.1 and 15% of ""^Ni positrons can 
escape in deflagration or delayed detonation models of SNIa, 
as well as in core-collapse SNe models. In contrast, the es- 
cape fractions of "^^Ti positrons range from 30 to 100%. The 
positrons that survived then enter the ISM where their slowing- 
down will proceed on longer timescales. The typical lifetime of 
MeV positrons in the ISM is about lO'^-lO'^yr, depending on 
the ISM phase (Jean et al. 2006 1, which allows their popula- 
tion to build up from many successive SNe. From their results 
on ''''Ni and "*"*Ti, suppler nented by an obser vational estimate of 



the positrons from ^^Al, Chan & Lingenfelter] (|l993jl showed 



that a steady-state production (and annihilation) rate of a few 
10*^ e^ s"' could be achieved, in agreement with the (model- 
dependent) value derived from observations. The study also il- 
lustrated the critical role played by the magnetic field topology 
inside the ejecta (which could confine positrons in the ejecta if it 
is strong and turbulent enough, thereby favouring energy losses 
and annihilation), and to a lesser extent by ejecta mixing (which 
could lift iron to the outermost lowest density regions of the 
ejecta, thereby favouring positron escape) and explosion energy 
(which controls the rate at which the ejecta becomes transparent 
to positrons). For instance, with the standard W7 deflagration 
model of |Nomoto et al.| ( |1984[ l, the authors found that the es- 
cape fraction of ^"Ni positrons in a SNIa ranges from 0.1 to 5% 
depending on whether the magnetic field is thorougly tangled 
or fully combed out, and this increases to 2.5 to 13% when the 
ejecta is uniformly mixed. 

Another way of estimating the escape fraction of decay 
positrons is based on the study of SNIa late lightcurves. The 
bolometric lightcurve of a SNIa is indeed powered by ^^Ni 
positron energy deposition after the ejecta became optically thin 
to the ^''Ni gamma decay photons, typically for t > 200 days. 
The driving parameters of the problem are again the structure 
and strength of the magnetic field inside the ejecta but also 
the ionisation state of the ejecta because an ionised medium 
is more efficient at slowing charged particles than is a neutral 
medium. Depending on whether positrons are trapped locally in 
the ejecta by a strong turbulent magnetic field or whether they 
can flow out almost freely along a weak radial magnetic field, 
the simulated bolometric lightcurves are different. In particular. 



deviations from one or the other scenario are detectable at late 
times, between 400 and 1000 days after explosion. From a set 
of models, Milne et al. ( 1999) derived positron escape fractions 
in the range ~0 to 11%, in agreement with the values found by 



Chan & Lingenfelter 



1993 1 (who did not include ionisation ef- 



fects in their study). From a comparison with a sample of 10 
lightcurves from a variety of SNe la (derived mostly from B- 
and V-band photometry), the authors concluded that positron es- 
cape occurs in SNIa ejecta and that a radial or weak magnetic 
field configuration is indicated by the observations. For the stan- 
dard W7 deflagration model and the assumption of a radial mag- 
netic field, they computed an escape fraction in the range 1.8 to 
5.5%, depending on the ionisation stage of the ejecta, with the 



higher value being favoured by the observations. Yet, Lair et al. 



(2006 1 showed from BVRI photometry that current data on late 
lightcurves of SNIa are insufficient to assess if color evolution 
and shift of the luminosity to longer wavelengths occur in the 
SNe la. Consequently, late bolometric light-curves extrapolated 
from optical bands could be increasingly inaccurate with time 
and would thus jeopardise the evaluation of the positron escape 
fraction. 

A third way to determine the escape fraction was explored 



by Kalemci et al. ( 2006) and is based on the annihilation emis- 
sion from SN1006, the remnant of a SNIa. The authors first used 
the INTEGRAL/SPI data to obtain an upper limit on the 5 1 IkeV 
flux from this object. Then, under the assumptions that escaped 
positrons can be locally confined close to the SNR and that their 
mean lifetime is less than 10^ yrs, they derive an upper limit on 
the escape fraction of 7.5%. These two strong assumptions on 
the diffusion and lifetime of positrons in the ISM (or CSM) are, 
however, not supported by any observational or even theoreti- 
cal argument. We will thereafter present a more detailed treat- 
ment of the annihilation of ^^Ni positrons in the ISM (and in the 
ejecta), which will allow a more accurate interpretation of the 
upper-limit on the 511 keV flux from SN1006. 

Most of the above-mentioned studies of positrons from SNe 
and SNRs were mainly focused on evaluating if these particles 
can escape and then fill the ISM and eventually feed the galactic 
annihilation emission. Yet, the direct annihilation emission from 
SNe and SNRs has been poorly studied. The simulations of SNIa 
late lightcurves including the contribution from decay positrons 
could have been the opportunity to provide simulated annihila- 
tion lightcurves from SNe la and their early remnants, but the 
limited time range of these simulations (up to ~ 1000 days) and 
the absence of constraints from gamma-ray line observations 
probably prevented the effort to be continued in that direction. 
INTEGRAL/SPI opened the way for cartography of the anni- 
hilation emission of diffuse or point-like origin. In KnodlsederJ 



et al. ( 2005 ), the authors exploited these capabilities to analyse 



the 511 keV emission from the Galaxy, including its extended 
component and the possible contributions from a multitude of 
point-like objects like SNRs. By now, the volume of available 
SPI data is as much as 6 times larger than the data set used by 
[Knodlseder et al.|P005| . Conditions are therefore favourable for 
a new search for annihilation emission from SNRs. 



4. INTEGRAL/SPI observations 

4.1. Instrument characteristics 

The SPI instrument onboard the INTEGRAL gamma-ray space 
observatory is a high-resolution spectrometer with imaging ca- 
pabilities (Vedrenne et al. 2003; Roques et al. 2003). The spec- 
trometry of gamma radiations with energies between 20 ke V and 
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Table 1. Upper limits on the 511 keV flux in units of 10 ^ phcm ^ s in three different energy bands and for various source sizes 
for SN1006 and Tycho. Also listed are the INTEGRAL/SPI exposure for each object in cm^ s. 





SPI exposure 


Sky model 


sin 1 9 Vf^\i 


JUo.J-JI J.J Kc V 
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point- source 
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3° gaussian 
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point- source 
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point- source 
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3° gaussian 
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6° gaussian 
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9° gaussian 
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12° gaussian 




7.9 








15° gaussian 
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SN1987A 


~ 8 X 10^ 


point- source 


3.8 


5.4 


6.6 


Gl.9+0.3 


~ 9 X 10^ 


point- source 


1.8 


2.6 


3.1 



8 MeV is performed by an array of 19 high-purity germanium 
detectors (as yet 17, since detector 2 and 17 failed on revolution 
142 and 210 respectively), with a spectral resolution of about 
2 keV FWHM at 511keV. All gamma-ray sources, diffuse or 
point-like, emitting in this energy range can also be imaged indi- 
rectly thanks to a coded mask system with an angular resolution 
of 2.8° and a fully-coded field of view of 16°xl6°. 

The data produced by SPI are strongly background- 
dominated (signal-to-noise ratio of less than 1 %) because sub- 
stantial hadronic interactions and electromagnetic cascades are 
induced by solar and cosmic-ray particles impacting the satel- 
lite. These result in an instrumental background composed of a 
continuum and a certain number of lines of various intensities. 
In particular, the 511 keV line appears in the background spec- 
trum and is even one of the strongest lines. Extracting the weak 
gamma-ray line signals of astrophysical origin therefore requires 
first an accurate modelling of that background noise. Then, the 
reconstruction of a given source intensity distribution is achieved 
through multiple exposures (following a dithering pattern) and 
model-fitting approaches (as explained below). 

4.2. Data set 

The data used for this analysis are all-sky data that have been col- 
lected between revolution 7 and 764 of the INTEGRAL satellite 
(the duration of one INTEGRAL orbit is ~3 days). A filtering 
has been applied to exclude all pointings during which abnor- 
mally high count-rates (due to solar flares and periodic radiation 
belt crossings) occurred and consequently preclude the detection 
of celestial signals. The total effective observation time eventu- 
ally amounts to 91.8 Ms, among which 58.6 Ms cover the galac- 
tic plane (|b| < 20°) and 33.2 Ms are high-latitude pointings (|b| 
> 20°). All single events dat^with energies in the 506-5 16 keV 
range have been selected and binned into 0.5 keV-wide bins. 



4.3. Method 

Due to the complexity of the response function, the SPI data 
cannot be inverted and the extraction of astrophysical signals is 
done by model-fitting. This implies a model for the instrumental 
background noise and a model for the sky intensity distributions, 
each model being composed of several component in the general 
case. 

In almost all gamma-ray line analyses performed from SPI 
data so far, the modelling of the background noise is achieved 
through combinations of activity tracers that are hypothesised 
to reproduce the various trends and time evolutions of the back- 
ground (for more details, see |Jean eTa l. 2003 ; Knodlseder et al. 
I005l [Martin et al.||2009b| 



a I. These tracers are strongly corre- 



lated with the data, but since we are looking for signals that 
are less than 1% of the data, they need to be finely adjusted 
through fitting. This approach has proven successful f or the 
study of the nuclear decay Unes from ^^Al, ^°Fe and "^"^Ti (Wang 
et al."2007] [Martin et aL]|2009b|a| , as weU as for the study of 
the diffuse 511 keV galactic emission ( |Knodlseder et al. [2005 [ 
Weidenspointner et aLp OOST). In the present work, we took ad- 
vantage of these past achievements and we therefore used the 



same background modelling as that employed by [Knodlseder 
([2^5^ and Weidenspointner et aL] ( [2008a| l 



etal. 



The sky models are assumptions on the intensity distribu- 
tion of the signal under investigation. They can be any kind of 
observed distribution that the searched emission is expected to 
be correlated with, or simple analytical models that are astro- 
physically meaningful and compatible with the performances of 
the instrument (on angular resolution especially). Each of these 
intensity distributions is mapped into the data-space through a 
convolution with the instrument response function (IRE) for each 
pointing of the data set and then forms a sky model component. 



Interactions of gamma photons with the detectors can be of two 
types: single events (SE), when the incoming photon deposits energy in 



one detector only, and multiple events (ME), when the incoming photon 
deposits energy in several adjacent detectors through Compton diffusion 
and pair-creation. 
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Table 2. Positron escape and survival fractions (in %) as a func- 
tion of parent nucleus, ejecta mass and transport scheme. 



Ejecta mass 


Transport 


56Ni 


■"Ti 


2 


Free 


0.4 


97.4 


2 


Trapped 


< 10-2 


90.8 


4 


Free 


< 10-2 


94.8 


4 


Trapped 


< 10-2 


79.7 


6 


Free 


< 10-2 


92.3 


6 


Trapped 


< 10-2 


69.0 


8 


Free 


< 10-2 


90.0 


8 


Trapped 


< 10-2 


59.2 



The sky model components and the background models are 
fitted simultaneously to the data using a Maximum Likelihood 
criterion for Poissonian statistics. The quality of the analysis is 
then assessed through examination of the residuals. Since we 
are working with background-dominated data, the modelling of 
the instrumental background can easily introduce systematic er- 
rors that might be compensated for by artificial signals from the 
sky. In addition, the assumption made on the sky intensity distri- 
bution might well be inappropriate or incomplete, which would 
introduce bias or leave unaccounted events in the data-space. 

To ensure that our models are adequate, we study the 
distribution of the residuals after back-projection on the sky. 
Basically, for all pointings of the data set, the residuals in each 
detector are uniformly projected across the mask (or more pre- 
cisely "across" the instrument response function) over the sky re- 
gion intersected by the field of view. The distribution of these sky 
residuals should then follow a statistical distribution, the shape 
of which is derived from simulated observations. This method is 
described in greater details in [Martin eFal . ( 2009b ) and proved 
for the present case that our hypotheses regarding both the back- 
ground noise and the sky intensity distribution are satisfactory 
and that the results obtained are not affected by systematic er- 
rors. 



4.4. Results 

In order to properly assess the 511 keV emission from the above- 
listed SNRs, we need to take into account the galactic diffuse 
emission at this energy to ensure that part of this diffuse flux 
is not attributed to the SNRs. This is particularly important for 
G1.9H-0.3, which lies very close to the galactic center, where 
most of the diffiise emission is found. The best-fit model so 
far for the diffuse annihilation emission is made of three com- 
ponents: a narrow gaussian distribution and a wide gaussian 
distribution (with typical angular sizes of 3° and 11° respec- 
tively and both centered on the galactic center) and a thick disk 
( Weidenspointner et al.^2008b) .We thus fit our models for the 6 
SNRs simultaneously to the 3-component model for the diffuse 
emission. 

Given the ~ 3° angular resolution of SPI, all 6 SNRs are mere 
point-sources for the instrument and we indeed started our anal- 
ysis with 6 point-sources at the positions of the SNRs plus the 
3 components of the diffuse emission. Yet, since positrons could 
have escaped and propagated away since the explosion, we need 



to search for annihilation emission not only in the SNRs but also 
in their surroundings. The largest spatial scale to be considered 
is I = ct, where f is the age of the remnant and c the light-speed. 
For SN1006 and Tycho, this propagation effect could lead to a 
substantial increase of the extent of the annihilation emission: 
up to 15° for SN1006 and up to 6° for Tycho. We have therefore 
explored how the result is affected by the source size by replac- 
ing the point-source model for these two objects by 2D gaussian 
profiles of increasing widths. 

Last, we also explored the potential broadening of the 
511 keV line by analysing the flux in three energy bands of 
increasing size, starting with the ~ 2keV (FWHM) spectral 
resolution of the instrument at 511 keV as the minimum size. 
Significant line broadening could arise from Doppler effect due 
to the following reasons: 

1. Annihilation in an expanding medium like the expanding 
steflar ejecta of a SN/SNR; for a velocity of ~ 1000 km s-', 
the line broadening is ~ 3 keV. 

2. Annihilation in a hot medium like the shocked regions of 
SNRs; for a temperature of 10^ K, the broadening is ~ 
30keV. 

3. Annihilation through formation of positronium in flight, 
which is an important process at the galactic scale and is 
expected to dominate in the neutral ejecta of SNe; the resul- 
tant broadening is ~ 6 keV ([Guessoum et al.|2005j|Jean et al. 



|2006| l. 

The spectrometric performances of the SPI instrument allow the 
study of the 511 ke V profile ( Jean et al.|2006^ but the corollary 
of such a capability is that too broadened signals become unde- 
tectable. 

The various sky models presented above have been fitted to 
the allsky SPI data, simultaneously to the background model. 
No 511 keV signal was detected from any of the 6 SNRs, what- 
ever the energy band or the spatial extent of the source, in agree- 
ment with the latest 511 keV cartography. The 3-component sky 
model for the diffuse emission proved robust and was not altered 
by the addition of the 5 SNRs, especially Gl. 9-1-0.3 which could 
have impacted the output due to its proximity to the galactic cen- 
ter. The added SNR sky models do not lead to any improvement 
of the fit and thus do not contribute to give a better representation 
of the SPI data. Finally, we note that using an asymetric model 
for the galactic diffuse emission (where the disk is split longitu- 
dinally in four parts, see [Weidenspointner et al. 2008b) gives the 
same results and does not modify the asymetry observed in the 
inner Galaxy. 

From our results, we derived 2cr upper-limits on the 511 keV 
flux from the 6 SNRs, for the various source sizes and line broad- 
ening assumptions. These are listed in Table [T] The values re- 
flect the uneven exposure of the sky by INTEGRAL. SN1987A 
is clearly the least exposed object among the 6 SNRs while 
Gl. 9-1-0. 3 is the most exposed. As expected, the limits become 
less constraining with increasing energy band or source size 
(because the signal is then diluted in a larger volume of the 
background-dominated data-space). 



5. Simulations 

To interpret the above observational constraints, we simulated 
annihilation lightcurves for a simple description of SN/SNR evo- 
lution and a simple uniform CSM/ISM, and under various pre- 
scriptions for the transport of the decay positrons. Many analyt- 
ical formulae developed here were already presented in [Chan &[ 
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|Lingenfelter| ( | 1 993) , hereafter CL93, and were just adapted for 
the purpose of computing annihilation Ughtcurves. 

5.1. SN/SNR model description 

In our model, the SN/SNR is represented as a uniform spherical 
ejecta in homologous expansion. In such a case, the velocity of 
the stellar material is a linear function of the radius and reaches 
a maximum v^j at the edge. The evolution of the ejecta is then 
given by the following relations: 



Rej - Vejt - 



3 M, 



10 E,j 
3 M,j 

3 



1/2 



Pej 



Anivejty 47r\10 



3 \3/2 M 
■J \ f 



5/2 



(2) 
(3) 




80 30 
Time Cyps) 



where M^j and E^j are the ejecta mass and kinetic energy, which 
are parameters of the model. R^j and p^j are then the ejecta size 
and density as a function of time t. 

The decay positrons of ^''Ni and '*'*Ti are injected at the very 
centre of the ejecta, where products of explosive Si-burning such 
as ^*Ni and ^"^Ti are expected to be found (yet, strong turbulence 
during the supernova can alter the ordered chemical stratification 
that results from a purely spherical explosion). We then consider 
two extreme transport schemes, as in 'CL93I either the positron 
is trapped in the ejecta, as a result of a strong magnetic turbu- 
lence for instance, or it is completely free and can travel radially 
up to the surface of the ejecta and then escape, provided it is not 
stopped before. 

Positrons lose energy along their way through ionisa- 
tion/excitation of the stellar material assumed to be neutral (but 
see|5.4|i. The energy loss-rate is given by the following law: 



dE _ Anrlm,c^p /Z 
dt u/3 



I I E + meC^X 



(4) 



where E is the kinetic energy of the positron, mgC~ its mass en- 
ergy and y - {\ -yS^)"'''^ its Lorentz factor. The characteristics 
of the medium are the density p, the mean ratio of atomic num- 
ber to mass number (Z/A) and the mean ionisation potential (/). 
The constants ro and u are the classical electron radius and the 
atomic mass unit. The term (l>{E) depends on the particle energy 
only and its complete expression is given in CL93| In our model. 



the medium is assumed to be of a homogeneous chemical com- 
position, with Z/A = 0.5 and (/) - 100 eV, both values being 
characteristic of the a-nuclei that make up most of the ejecta of 
SNe. 

If we rewrite Eq.|4]in a more general wa}|^ 



dE 
dt 



= -p<D(£) 



pdt — - 



dE 
0(£) 



(5) 



then the initial and final states of a slowing-down phase, respec- 
tively (f,, Ei) and (f/, £/), can be related through: 



^! dE 



f(£) 



(6) 



The integrated density history experienced by the particle be- 
tween f, and tf is thus related to an integral in the energy space 



' The function Q>{E) actually has a dependence on the chemical com- 
position of the medium (through Z/A and I) but since we work with a 
homogeneous chemical composition, we dropped that dependence. 
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Fig.l. SllkeV lightcurve of a SN/SNR for and 
positrons (upper and lower panel respectively) annihilating in 
the ejecta, for various ejecta masses and the two transport sce- 
narios. 



between E-, and Ef. Combining Eq. |6|with Eq. [3j it is possible 
to compute the time necessary for aireshly injected particle to 
be slowed down to a certain threshold energy E,hr, taken in the 
following to be 100 eV: 



K 
2 



(I 



Je,i, 



dE 
0(£) 



where p^y - Kt 



(7) 



In the continuously decreasing density of the expanding ejecta, 
there is for each injection time a maximum energy above which 
the particle cannot be slowed down to E,hr. This energy limit 
£■;;,„ decreases as time passes and is obtained by taking tsd — > °° 
in the above equation: 



K 1 _ r^'" 



dE 



(8) 



Practically, this means that an increasing fraction of the positrons 
released by the decay of ^''Ni or '^"^Ti will not be slowed down 
enough and hence will not annihilate in our model. This is due 
to our assumption of homologous expansion for the SN/SNR 
evolution, which allows the density to go to zero at very late 
times. In reality, however, homologous expansion is halted by 
the inward propagation of a reverse-shock that basically prevents 
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the ejecta density from falling far below that of the surrounding 
CSM/ISM. In addition, positrons may escape the SNR and en- 
ter this CSM/ISM where their slowing-down proceeds within a 
finite timescale. 

In the context of decay positrons in SNe/SNRs, the initial 
state (f,, £/) of a positron is described statistically: the injec- 
tion time follows an exponential distribution characterised by the 
half-life of the parent nucleus, while the initial energy obeys a y6- 
spectrum that depends on the charge of the parent nucleus and 
on the endpoint energy of the positron (the mathematical form 
of the ;0-spectrum is given in CL93[ l. For each initial state, we 
compute the slowing-down time tsd using Eq.[7] together with an 



estimatqjof the crossing time tc-oss = vti/{v - Vej) at which the 
particle reaches the surface of the ejecta (with v being the particle 
velocity). Another relevant time is f j„„„, which is the range over 
which the annihilation lightcurves are computed and is set here 
at llOOyrs. Depending on the obtained values and the adopted 
transport scheme: 

1. Trapped positrons: If f ,^ < f5,m„, the particle is slowed down 
within the typical lifetime of a young SNR and it annihilates 
at a time tsd- 

2. Trapped positrons: If tsd > tsimu, the particle cannot be 
slowed down and it survives inside the ejecta. 

3. Free positrons: If t^d < t cross ^ tsi„„„ the particle is slowed 
down before it can escape the ejecta and it annihilates at a 
time tsd- 

4. Free positrons: If tsd > tcross, the particle can escape into the 
CSM/ISM the ejecta before it is completely slowed down. 

In the above scenario, we considered that once slowed down to 
the threshold energy, the positrons all annihilate through positro- 
nium formation in-flight in the ejecta assumed to be mostly neu- 



tral (see 2.2 1. The process of positronium formation and sub- 
sequent annihilation is considered to be immediate owing to its 
very large cross-section and the ns-/is lifetime of the bound-state. 

The model described above yields the time distributions 
of annihilation for ^^Ni or ''^Ti positrons, for a given ejecta 
mass and energy and for the two transport schemes. SN/SNR 
lightcurves at 5 1 1 ke V are eventually obtained after normalisa- 
tion by the following factor: 



fnorm ~ 



yAfpifps 1 1 

Au 4 AjkP- dt 



(9) 



where Ya and A are respectively the typical yield of the radionu- 
clei and its mass number, and fp is the branching fraction of 
yS^-decay for that isotope. Then, fps is the assumed positronium 
fraction and the factors 1/4 and 2 account for the fact that 25% of 
the positroniums are parapositroniums, which decay by emitting 
2 photons at 5 1 1 keV. Last, d is the distance to the SNR and dt 
is the size of time bins in our simulation. 

The transport of 5 1 1 ke V photons in the SN/SNR is taken 
into account in a very basic way, by comparing the mean free 
path of the photons in the uniform density ejecta to the radius of 
the ejecta. In the computation of the 511 keV photon mean free 
path, we used a value of 8.63 lO^^cm^g"' for the attenuation 
factor, which formally corresponds to the O element but also ap- 



* The formula used to estimate f^.„.„ assumes that the particle velocity 
remains constant at its initial value. While moving through the ejecta, 
the particle velocity decreases so the true crossing time is larger than 
what we compute. Yet, as shown in |5.3| most escaped particles are still 
relativistic when they exit the ejecta so the error on the crossing time 
can be considered as small. 




0,4 0,6 0,8 1 

Energy (MeV) 



Fig. 2. Spectra of the escaping ^''Ni and '^'^Ti positrons (upper 
and lower panel respectively), as a function of ejecta mass. Also 
shown for comparison is the original y6-spectrum. 



plies to elements from Ne to n]^ 

In the process, we also computed the spectra of the escaped 
particles as these will be the input to the modelling of the an- 
nihilation in the surrounding CSM/ISM. The energy E^sca of an 
escaped particle is given by: 



(\ 



dE 



(10) 



As we will see later, escaped particles in our model have a mean 
kinetic energy that is a few 100 keV below the mean kinetic en- 
ergy at decay, which contributes to an enhanced annihilation in 
the suiTounding CSM/ISM. 

5.2. CSM/ISM model description 

The model for the CSM/ISM is actually a simplified version of 
that used for SNe/SNRs. The density in the medium surrounding 
the SN/SNR is assumed to be constant, which results in a more 



The attenuation factor values were taken from the Xcom 
database of the National Institute of Standards and Technology; see 
"http://physics.nist.gov/PhysRefData/Xcom/html/xcoml.htmr' 
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simple equation to compute the slowing-down time: 



Po(tsd - ti) 



dE 



(11) 



Then, the methodology to compute annihilation lightcurves is 
the same as for the case of SNe/SNRs, except that an escape 
fraction is added to the normalisation factor of Eq. |9] Another 
difference is that the energy distribution of the positrons injected 
in the CSM/ISM may differ from the j6-spectrum that prevailed 
inside the SN/SNR because the particles already lost a fraction 
of their initial energy on their way out. 

If the CSM/ISM is assumed to be neutral, energy loss pro- 
ceeds through ionisation/excitation and annihilation again oc- 
curs mainly by positronium formation in-flight immediately af- 
ter slowing-down. If the CSM/ISM is assumed to be strongly 
ionised, energy loss proceed through Coulomb losses, which are 
~10 times more efficient than excitation/ionisation and would 
therefore favour earlier annihilation. Yet, if the CSM/ISM is 
assumed to be fully ionised, no H atoms exist for positron- 
ium formation in-flight and annihilation occurs mostly through 
positronium formation by radiative recombination of ther- 
malised positrons with free electrons. This process, however, 
has quite a small cross-section and takes place over character- 
istic timescales of ~ 10^ n^' yrs (where is the electron den- 
sity), which strongly limits the annihilation in the vicinity of 
the SNR and rather favours diffusion out to large distances. In 
a warm ionised medium, cosmic interstellar He is essentially 
neutral (Ferriere 2001] ), but positronium formation by charge 
exchange with He is negligible as positrons are preferentially 
slowed down below the energy threshold of the process by the 
free electrons. Similarly, annihilation in dust grains and on PAHs 
contributes only marginally mainly because of their very low 
densities ( |Guessoum et al. |2006| ). On the whole, it seems that 
a strongly ionised medium is not a favourable condition for 
annihilation on timescales of the order of the typical ages of 
young SNRs, and we will therefore consider in the following 
only slowing-down and annihilation in a predominantly neutral 
CSM/[SM. 



5.3. Results 

Based on the above model, we computed SN/SNR lightcurves at 
511 keV for both ^''Ni and '^'^Ti positrons, for a single canonical 
explosion energy of 10^' erg and a range of ejecta masses, and 
under the two extreme prescriptions for positron transport. In or- 
der to allow an easy scaling of our results (so as to apply it to 
specific cases), we used isotope yields of 10"' M© and 10^* Mq 
for ^^Ni and '*'*Ti respectively, together with an assumed positro- 
nium fraction of 90% and a distance to the SN/SNR of 1 kpc. 
The results are presented in Fig.[T] 

Whatever the transport scenario, the annihilation lightcurves 
for ^^Ni positrons quickly rise (within a year) to quite high 
fluxes, with a peak value of ~l-10phcm"^ s ' (for a distance 
of 1 kpc), that is about 10^-10^ times the present-day sensitivity 
of INTEGRAL/SPI for point-sources (see Table [T]i. After that 
"annihilation flash", the emission rapidly drops. Over the first 
years, the lightcurves closely follow the decay of ^^Ni, which 
can be explained by a short slowing-down time with respect to 
the characteristic time of the decay. Then, the lightcurves depart 
from the decay curve as the lifetime of the particles becomes 
long enough for an increasing fraction of the free positrons to 
escape, while the trapped positrons annihilate with an increas- 





Fig.3. 511 keV lightcurve of a SN/SNR for ^^Ni and ^^Ti 
positrons (upper and lower panel respectively) annihilating in 
the surrounding CSM/ISM after escape, for various atomic hy- 
drogen densities. The energy distribution of the particles is their 
original yS-spectrum. 



ing delajj^ Lighter ejecta favour an early escape and the cor- 
responding annihilation lightcurves for free positrons therefore 
drop more rapidly. In the case of trapped positrons, the lower 
densities associated with the very lightest ejecta imply longer 
slowing-down times and hence a higher flux at later times. 

The annihilation lightcurves for "^"'Ti seem to follow similar 
trends. Nevertheless, the continous injection of ^'^Ti positrons 
over centuries, as opposed to the burst of ^^Ni positrons, leads to 
some differences. In particular, the delayed emission of trapped 
positrons increases with ejecta mass, contrary to the case of ^^Ni 
positrons. In the latter situation, almost all positrons are released 
at early times and need to survive as long as possible for substan- 
tial late emission to exist, which requires light and thin ejecta; in 
contrast, most '*"'Ti positrons are released at intermediate times, 
when the ejecta is already considerably diluted, and the key for 
late emission is then an efficient slowing-down, which requires 



In the trapped case, one may have the impression that the 
lightcurves are non-physical since after following the decay curve, they 
rise above it thereby suggesting that more positrons are annihilated 
than created. This is a misleading effect of the many orders of mag- 
nitude spanned by the lightcurves. Actually, only a tiny fraction of the 
positrons created over the first years annihilate with a delay and the 
delayed emission a few decades after the explosion is accordingly tiny 
compared to the flux levels of the early years. 
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Table 3. Characteristics of the six SNRs used to compute their 511 keV annihilation lightcurves. Numbers in italic are typical values 
assumed when no established measurements were found in the literature. 





Cas A 


Tycho 


Kepler 


SN1006 


SN1987A 


Gl.9+0.3 


Age (yrs) 


338 


437 


405 


1003 


22 


100" 


Distance (kpc) 


3.4'' 


2.4'^ 


6" 


2.2 = 


50 


8.5" 


Type 


SNllb f 


SNla 


SNla 


SNla 


SNll 


SNla" 


Ejecta mass (M©) 


2.2- 


1.4 


1.4 


1.4 


14 h 


1.4 


Explosion energy (10" erg) 


2.0-4.0 ' 


1 


1 


1 


1.1 ±0.3 


1 


^''Ni yield (Mq) 


0.07-0.15 f 


0.6 


0.6 


0.6 


0.07 


0.6 


""•Ti yield (M^) 


1.6 lO-'* J 


< 2 lO-* 


< 3 10-5 ' 


< 3 lO-"' ' 


210-'' 


< 3 10-5 1 


CSM density (cm"^^) 


3.2" 


0.3™ 


0.5'' 


0.1 " 


0.02° 


0.04" 



Refer ences: (a) Reynolds et al. (20081 (b) Ree d'etari(fT995ll (c)|Cassam-Chenai et al.jpOOTl l (d) Vink (2008) (e) Winkler etari(|2003) (f)|Krause| 
eraLl([2008|l (g),Willingale et al.,(|2003|l (h),Blinrukov et al.H2000> (i) [Laming et al.|(|2006ir(j)[Renaud et al., ( |2006bj (k) |Renaud et al.| ( |2009^ (1) 
Renaud et air| ( |2006a^ (m) |Truelove & McKeeH1999^ (n) |Raymond et al.| ( |2007| ) (o) |Sugerman et al.| ( |2005^ 
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Fig. 4. 511 keV lightcurve of the CSM/ISM for ^^Ni or ^^Ti 
positrons from SNe/SNRs, for an atomic hydrogen density of 
1 cm-^ and depending on the energy distribution of the particles. 



thick and massive ejecta. On the whole, the most favourable con- 
ditions to observe 511 ke V emission from young Galactic SNRs 
therefore seem to be a massive ejecta, a high yield of '^"^Ti and 
enough magnetic turbulence so as to confine the decay positrons 
inside the ejecta. 

The escape fractions obtained in the case of free positrons 
are listed in Table |2] together with the survival fractions obtained 
in the case of trapped positrons. The ~ 0% escape fracti ons for 



""^Ni positrons agree with the negligible values found by |CL93 
for similar ejecta masses and no mixing (their models B4B, B6C 
and B8A can be approximately compared to our 2, 4 and 6Mq 
ejecta models). CL93 showed that higher escape fractions of a 
few % can be obtained if the ^^Ni is mixed out in the ejecta, for 
higher explosion energies, or in the case of lighter ejecta masses 
(of SNla for instance), but we did not explore these effects here. 

The spectra of the escaped free particles are shown in Fig.|2] 
for the two parent radio-isotopes (in the case of ^^Ni, the escape 
fractions are so small that the spectra are shown on a logarithmic 
scale, for the 2 and 4M0 ejecta only). The original y6-spectrum is 
also shown for comparison. In both cases, the number of escaped 



particles decreases as ejecta mass increases and their spectrum is 
shifted to lower and lower energies. This will have appreciable 
consequences for the annihilation in the CSM/ISM. 

The 511keV lightcurves for the annihilation of decay 
positrons in the CSM/ISM surrounding the SNR are shown in 
Fig. [3] for various atomic hydrogen densities. Again, in order 
to allow an easy scaling of our results, we used isotope yields 
of IQ-' Mq and IQ-"^ Mq for ^''Ni and '*''Ti respectively, together 
with an assumed positronium fraction of 90% and a distance to 
the SN/SNR of 1 kpc. In addition, the escape fractions are as- 
sumed to be 1% for ^''Ni positrons and 100% for ^"^Ti positrons. 

In Fig. [3] the positrons were injected with an energy distri- 
bution corresponding to their original and unaltered jS-spectrum 
and the lightcurves obtained for ^''Ni and "''*Ti positrons are re- 
markably similar The fact that nearly identical levels of emis- 
sion are reached can actually be explained by the mean charac- 
teristics of each positron source: although the yield of ^^Ni is 
on average 10-* higher than the yield of '''^Ti, only 20% of the 
''^Ni nuclei undergo a yS^-decay and because of the short life- 
time of the isotope only ~ 1% of the positrons can escape; in 
contrast, almost all '*'*Ti nuclei produce positrons and almost all 
these positrons escape. 

It should also be noted that the evolution of the annihilation 
lightcurves with density agrees with expectations. If a flux /511 
is obtained at t for a CSM/ISM density «i, a flux («2/«i)/5ii 
is obtained at {ni/n2)t for a CSM/ISM density n2 (because the 
slowing-down time decreases as n increases, and so successive 
positron annihilations in a higher density medium occurs sooner 
and in a smaller time interval). The scaling works as expected for 
^^Ni, but not so well for '''^Ti due to its longer decay time (which 
prevents, for instance, the comparison of the flux at 1000 yrs in 
the « = 0.1 cm--' case with the flux at 100 yrs in the n = LOcm--* 
case because all positrons have been released at 1000 yrs but not 
at 100 yrs). 

Prior to entering the CSM/ISM, however, most positrons 
lost a fraction of their initial energy on their way through the 
ejecta (see Fig. |2]). When these modified spectra are taken into 
account, the resulting annihilation in the CSM/ISM is notice- 
ably different. In Fig. [4] the lightcurves obtained for ""^Ni and 
"'"'Ti positrons entering a 1 cm--' CSM/ISM after crossing a 4 M© 
ejecta are shown, together with the unaltered yS-spectrum case (a 
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single curve is shown for both isotopes since, as discussed ear- 
lier, they are quite similar). It appears that the annihilation flux 
from the CSM/ISM can be increased by up to two orders of mag- 
nitude in the case of ^^Ni positrons when the slowing-down of 
the particles inside a 4Mq ejecta is taken into account. We want 
to emphasise, however, some simplifications made when com- 
puting the CSM/ISM lightcurves with modified positron spec- 
tra. First, we used escape fractions of 1% and 100% for ^''Ni and 
^'^Ti positrons respectively, but these are not the values given by 
our model for a 4Mo ejecta (see Table |2]i. We did that to ease 
the comparison with Fig. |3] but one should keep in mind that the 
escape fraction and the spectrum of escaped particles are linked 
(higher escape fractions are associated with more energetic es- 
caping particles). Then, we used the same modified spectrum 
for all positrons injected in the CSM/ISM, but this spectrum ac- 
tually evolves with time (the mean energy of escaping particles 
increases with time, as the ejecta becomes more and more trans- 
parent). Nevertheless, the comparison of simulations made for 
different modified spectra and different injection time profiles 
showed that these simplifications only have a modest impact on 
the predicted lightcurves. 
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Fig. 5. Adiabatic losses and ionisation/excitation losses as a 
function of time for a 2 Mq ejecta and for 1 keV and 1 MeV par- 
ticle energies. 



5.4. Adiabatic losses 

In our model of positron transport in expanding SN/SNR ejecta, 
we assumed that the particles lose energy through ionisa- 
tion/excitation of the stellar material. Theoretically, however, 
positrons also suffer from adiabatic losses. In the following, we 
discuss how these additional losses may impact the predicted 
lightcurves. 

The general expression for adiabatic losses is given by Eq. 
12 where E is the energy of the gas particles, y the ratio of spe- 
cific heats and v the velocity of the flow. 



dE 
dt 



= -(r-l)(V.v)£ 



(12) 



In the specific case of homologous expansion (v - rj t) of an 
ultrarelativistic gas (y - 4/3), this equation becomes: 



dE 
dt 



E 
t 



(13) 



In contrast to energy losses by ionisation/excitation, energy 
losses by adiabatic expansion do not depend on the character- 
istics of the ejecta (such as mass or kinetic energy) and they 
increase with particle energy. In order to be efficiently slowed 
down by adiabatic losses, however, the particles need to be con- 
fined to the expanding volume for a sufficiently long time com- 
pared to the expansion timescale. 

In Fig.|5]are plotted the energy losses by both processes as a 
function of time, for an ejecta mass of 2 M© and for particle en- 
ergies of 1 keV and 1 MeV. From this comparison, it appears that 
adiabatic losses can be neglected over the first years after explo- 
sion. Then, owing to the fact that ionisation/excitation losses de- 
crease as while adiabatic losses decrease as 1/f only, the lat- 
ter become progressively dominant. Energy losses by adiabatic 
expansion take over after ~5-10 yrs for ~MeV particles and after 
~500-1000yrs for ~keV particles. For higher ejecta masses, the 
shift occurs later as the associated higher densities (for a same 
explosion energy) maintain ionisation/excitation losses over adi- 
abatic losses for a longer time (ionisation/excitation losses scale 
with density, which scales as the power 5/2 of ejecta mass; see 
Eq.[3]). 

In the context of our simple model, it is not possible to 



include both processes in Eq. |6] to compute the slowing-down 
time because of their different time dependence. Nevertheless, 
we show below from simple arguments that adiabatic losses are 
not expected to modify the annihilation lightcurves computed 
from our model. 



In 5.1 we showed that there exists for each particle injected 
in the ejecta at time t an energy limit Eumit) above which the 
particle cannot be slowed down in a finite time (see Eq. |8]l. 
Under the assumption of an undisturbed ejecta expansion, this 
energy limit gives the increasing fraction of trapped positrons 
that will survive in the ejecta, mostly as relativistic particles. 
It also gives to a very good approximation the increasing frac- 
tion of free positrons that will escape the ejecta. Although not 
really intuitive, it can be easily shown that for each time t af- 
ter explosion, the particle energy above which adiabatic losses 
dominate over ionisation/excitation losses has almost the same 
expression as £/,„, (one just need to equate both energy loss for- 
mulae to get the result). In a very opportune way, adiabatic losses 
become important at a given energy when ionisation/excitation 
losses fail to efficiently slow down positrons with that initial en- 
ergy (and this works whatever the ejecta mass and explosion en- 
ergy). At this point, the energy loss rate of adiabatic expansion 
and ionisation/excitation are almost equal, but because ionisa- 
tion/excitation losses decline faster than adiabatic losses, the par- 
ticle rapidly suffers from adiabatic losses only. From integration 
of Eq. [13] the slowing-down time for a particle injected at time 
t, with initial energy E, and suffering adiabatic losses is found to 
b^: 

£■ 

tsd = TT^ ti (14) 

From Eq. [14] and Fig. [5] we can now estimate how positrons 
suffering adiabatic losses can affect the annihilation lightcurves 
for the case of a 2Mq ejecta. Positrons with initial energy of 
~ 1 MeV (respectively ~ 1 keV) are affected by adiabatic losses 
from 5-6 yrs (respectively 500-600 yrs) after explosion, and they 
are slowed down to £■,;„• after 50000-60000 yrs (respectively 
5000-6000 yrs). Adiabatic losses in a 2Mq ejecta expanding 



Strictly speaking, Eq. 



14 



should include a numerical factor of the 
order of a few because the energy loss rate of Eq. [13] increases by up 
to a factor of 2 as particles turn from ultra-relativistic to mildly or non- 
relativistic. 
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with a kinetic energy of 10^' erg therefore lead to annihilation 
at very late times, well beyond the typical ages of the SNRs in 
our sample. Moreover, light ejecta are the most favourable case 
because adiabatic losses take over more rapidly, and the corre- 
sponding annihilation occurs at the soonest. We can thus safely 
conclude that adiabatic losses would impact our predicted anni- 
hilation lightcurves only marginally. 



6. Discussion 

We compared our upper-limits on the 511 keV emission from 
the youngest local SNRs with the expected fluxes computed in 
the frame of our simple models. In this comparison, we took 
into account the peculiarities of each object that can impact the 
annihilation emission. 



6. 1 . Characteristics of the young SNRs 

The main parameters required to model the annihilation in an 
individual SNR are the ejecta mass, the ^^Ni and '^'^Ti yields, the 
density of the surrounding medium and the distance to the ob- 
ject. In Table[3]are listed the relevant characteristics of the SNRs 
considered in this study. The quantities for which we did not find 
clearly established measurements in the literature were assumed 
to have typical values. 

The ejecta mass for SNIa explosions was set to 1.4 Mq, 
with the implicit assumption that these events correspond to 
the deflagration or detonation of a white dwarf that reached the 
Chandrasekhar mass through accretion of matter from a compan- 
ion star. For Cas A, the ejecta mass was determined from mod- 
elling of the X-ray emission of the remnant, while for SN1987A 
it was determined from modelling of the optical/UV lightcurve 
of the supernova over the first months. 

The -''^Ni yield of afl SNe la was set to 0.6 Mq, as deduced 
from the quite homogeneous early lightcurves of this type of 
events. The ^^Ni yield of Cas A was assumed to be in the range 
inferred for SN1993J because the optical spectrum near peak 
brightness of Cas A, echoed by interstellar dust several centuries 
after the explosion, is quite similar to the optical spectrum of 
SN1993J, which is taken to be a canonical SNllb. In contrast, 
the ^''Ni yield of SN1987A was deduced dkectly from its early 
lightcurve. 

The production of ^"^Ti by SNe remains quite controversial 
because Cas A is currently the only observed source of gamma- 
ray radiation from "*'*Ti decay, while emission from younger and 
probably obscured SNRs is statistically expected (see the discus- 
sion in The et al.,2006 ). The '''*Ti yield of the Cas A explosion 
was estimated to 1.610~'*Mo from direct observations of the 
gamma-ray decay lines. A similar value was indirectly inferred 
for SN1987A from the modelling of the late lightcurve of the 
event. For Tycho and Gl. 9-1-0. 3, upper-limits on the flux at the 
energy of the two hard X-ray decay lines at 68 and 78 keV were 
translated into 3cr upper-limits on the mass of "^"^Ti synthesised 
in the explosions. These limits are consistent with most theoret- 
ical yields obtained in deflagration and delayed detonation mod- 
els of SNIa explosions, which range from 0.8 to 4.6 10"^ Mq 
( Iwamoto et al.| 1999 1. For all SNe la, we therefore adopted a 
yield of 3 10"^ Mq, in agreement with the predictions and the 
strongest constraint set by Gl.94-0.3. 

The CSM/ISM densities around SNe/SNRs can be estimated 
by various methods. From a modelling of its expansion dynam- 
ics, Cas A was found to be expanding in the dense wind of its red 



supergiant progenitor The same method was applied to Tychcp^ 
Kepler and Gl. 9+0.3. For SN1006, the 0.1 cm^^ CSM/ISM den- 
sity is imposed by the observed thickness of a Balmer-dominated 
filament of the forward shock. For SN1987A, the density we 
adopted corresponds to the average value obtained from the 
~2 Mq of material that were estimated from light-echoes to form 
the circumstellar environment of the progenitor within about 30 
light-years. 

Most of these works, however, rely on the assumption of a 
constant density environment while the actual density structure 
around a SNR is quite often inhomogeneous and evolves with 
distance. In the general case, the density structure around a core- 
collapse SN results from the interaction of successive stellar out- 
flow episodes, each with a specific mass loss rate and velocity. 
In the spherical case, this can give rise to a series of concen- 
tric shells of various densities and thicknesses ( |Garcia-Segura 
et al.||199 6a b ; Dwarkadas 2005 ), but reality seems to be even 
more complicated. Light-echoes from SN1987A have revealed 
the complex circumstellar environment shaped by the progeni- 
tor B-star before its explosion: the slow and dense wind blown 
during the red supergiant phase is thought to have been asym- 
metric and, while expanding in the matter expelled during Main- 
Sequence, created an equatorial overdensity that forced the sub- 
sequent blue supergiant wind to polar directions. The resulting 
double-lobed structure spans several orders of magnitude in den- 
sity and is the proposed explanation for the three rings that ap- 
peared shortly after the supernova (see Sugerman et al. 2005 and 
references therein). The case of thermonuclear SNe is not less 
complicated. During the pre-SN evolution, as the white dwarf 
grows in mass through accretion, the surrounding environment 
can be shaped by a variety of processes among which the mass- 
loss from a red giant companion, a fast accretion wind or recur- 
rent nova explosions ( Badenes et al.|2007 Hachisu et al.|2008 
Borkowski et al. 2009). 

In addition to the uncertainties on the density structure 
around SNe/SNRs, the physical state of the surrounding gas is 
also relevant to the annihilation of decay positrons. As explained 
previously, a strongly ionised medium is more efficient at slow- 
ing down the escaping positrons but then implies annihilation 
on very long timescales (compared to the typical ages of young 
SNRs). The massive star progenitors of core-collapse SNe are 
powerful sources of UV radiation that ionise the surrounding 
medium over large distances. In addition, the burst of soft X-rays 
and extreme UV at supernova shock breakout ionise or reionise 
the CSM/ISM up to certain distance. The timescale for recombi- 
nation in the CSM/ISM may be longer than the age of the SNRs 
and so escaping positrons may not annihilate efficiently over the 
first centuries. 



6.2. Observations versus predictions 

From the data listed in Table [3] we computed the expected fluxes 
for each SNR. As observational constraints, we considered the 
INTEGRAL/SPI upper limits on the 51 1 keV flux obtained in a 
5 keV band around the line, because annihilation in our model 
was taken to occur mostly through positronium formation in- 
flight, which causes a Doppler broadening of about 6 keV (see 
4.4 1. Moreover, we assumed the most constraining case of no 
spatial diffusion. In the following, we discuss separately the case 



We note that, for Tycho, the 0.3 cm ' density derived by 



Truelove 



|& McKee| l |199 9 1 from d ynamical arguments i s cons istent with the 



0.6 cm " limit obtained by 



Cassam-Chenai et al. 



(2007 



to account for 



the absence of thermal X-ray emission from the shocked ambient gas 
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of annihilation in the expanding ejecta and in the CSM/ISM. 

Instead of showing specific synthetic 511 keV lightcurves for 
each remnant of a SNla, we took advantage of the homogeneity 
of this class of object to model a single set of lightcurves for 
a paradigmatic SNla located at a distance of 1 kpc. Our upper- 
limits on the annihilation fluxes were then scaled to that distance 
and compared to the predictions. 

The resulting lightcurves, including the contributions from 
both ^^Ni and **Ti positrons and depending on the transport sce- 
nario, are shown in Fig. [6] together with the constraints obtained 
from INTEGRAL/SPI observations. The absence of 511keV 
emission from the remnants of recent Galactic SNe la turns out 
to be consistent with the predictions of our simple model. In the 
most favourable case of positrons being trapped in the ejecta, 
the emission starts at quite high values but falls by more than 7 
orders of magnitude over the first 50 yrs and then continues de- 
creasing more gently over the following centuries. Our observa- 
tions of Tycho, Kepler and SN1006 are definitely not constrain- 
ing (and are therefore not shown in the plot). The upper-limit on 
the 511 keV flux from the youngest Galactic SNR is 4 orders of 
magnitude above the predicted flux for the trapped scenario, and 
so Gl. 9-1-0.3 does not help constraining the process of positron 
transport and annihilation in SNe/SNRs either. 

For the two remnants of core-collapse SNe of our set, 
we computed specific annihilation lightcurves shown in Fig. [7] 
where the contributions from both ^''Ni and '^"'Ti positrons were 
added. In the case of Cas A, our upper-limit is more than 3 orders 
of magnitude above the highest predicted flux, obtained when 
positrons are trapped in the ejecta. The situation is similar for 
SN1987A, although our upper-limit is closer to the theoretical 
lightcurves. In both cases, the INTEGRAL/SPI observations are 
not constraining enough to favour one scenario of positron trans- 
port over the other 

As seen above, the highest annihilation fluxes from SNRs are 
obtained when the decay positrons are trapped in the expanding 
ejecta. When the latter are free to escape, the corresponding an- 
nihilation emission strongly drops within a few decades after the 
explosion. In this case, however, their slowing-down and annihi- 
lation proceed in the CSM/ISM. From our simulation of decay 
positron annihilation in the CSM/ISM, we found that escaped 
positrons can give rise to the following steady 511 ke V flux over 
the first millenary: 



-^511 -KejKe, 



1 kpc 



\lcm-3/ 



^56 
0.1 Me 



(15) 



X 10 ""phcm ^ s ' 



where d, no and Fsf, are respectively the distance to the object, 
the CSM/ISM density and the iron yield. K^sc is the escape frac- 
tion in % and K^j is a correcting factor accounting for the effect 
of positron energy loss due to ejecta crossing, which was found 
empirically to be about 1 for light ejecta of I-2M0 and about 2 
for massive ejecta of 8-1 M ^. The formula involves only ^''Ni 
because it was shown in 



5.3 



that ^"'Ti positrons escaping in the 
CSM/ISM lead to an annihilation flux that is more than a factor 
of 10 lower than the flux arising from escaping ^''Ni positrons 
(see Fig.|4]i. The contribution of "^'^Ti is therefore neglected. 

From Eq. [15] and the data in Table [3] we can compute the 
expected CSM/ISM annihilation fluxes and convert our upper- 
limits on the 511 keV flux into upper-limits on the escape frac- 
tions. We eventually obtain upper-limits on the escape fraction of 
-13% for Cas A, -12% for Tycho, -30% for Kepler and -33% 
for SN1006. SN1987A and Gl. 9+0.3 are too far away to be con- 
straining. In the case of substantial spatial diffusion of escaped 
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Fig. 6. Total 511 keV lightcurve from a SNla at a 1 kpc dis- 
tance, combining the contribution of ^^Ni and '^^Ti positrons. 
Also shown are the upper-limits on the 511 keV flux from 
INTEGRAL/SPI. 



positrons in the CSM/ISM, the upper-limits on the escape frac- 
tion in Tycho and SN1006 would be increased by a factor of -2 
(but in this case, our estimate for the surrounding density may 
not be valid anymore). 

The value obtained for SN1006 is higher than the 7.5% 
estimate from jKalemci et al.| ([2006 ) but was derived from a 
much finer model and should therefore be considered as more 
reliable. Our most constraining upper limits of 12 and 13% on 
the escape fraction of ^^Ni positrons in Tycho and Cas A, re- 
spectively, are consistent with the maximum values obtained 
from modelling by Chan & Lingenfelter ( 1993 ) and Milne et al 
(1999) and reported in Sect. [3] In addition, ,Higdon et al., ( ,2009| 
showed that the steady-state Galactic positron production rate 
of - 2.4 X lO'^^e+s ' and afl other properties of the 511keV 
INTEGRAL/SPI observations can be entirely explained by a nu- 
cleosynthesis origin if the ^*Ni positrons from SNe la have an 
escape fraction of 5 +2%. In this case, ^^Ni would provide about 
two thirds of the Galactic positrons, the remainder coming from 
''^Ti (20%) and -^Al (12%) for which escape from SNe/SNRs 
is not a problem, and these -MeV positrons can propagate over 
fairly large distances in the Galaxy before they annihilate. An 
escape fraction of 5% for ^^Ni positrons actually corresponds to 
the value obtained by [Chan & Lingenfelter| ( |1993| l for the W7 
deflagration model of Nom oto et al.| (1984) under the assump- 
tion of no ejecta mixing and a combed out magnetic field. It is 
also consistent with the 5.5% obtained by Milne et al. ( 1999 1 for 
the same W7 model, a radial or weak magnetic field, and a 1% 
ionisation of the ejecta. Our 12% upper-limit on the escape frac- 
tion of ^^Ni positrons for the Tycho SNla remnant is therefore 
consistent with these estimates, and therefore with a nucleosyn- 
thesis origin of the Galactic positrons that annihilate at 5 1 1 keV. 

We emphasise that the above results are associated with 
the assumptions of a uniform neutral ejecta homologously ex- 
panding in a uniform neutral medium. More realistically, super- 
nova ejecta rather follow exponential or power-law density dis- 
tributions (see f or instance [ Marcaide e t al.|1997) [Dwarkadas &| 
Chevalier' 1998] [Truelov e & McKee 1999) and their expansion 
is accompanied by the formation of shocks that heat and com- 
press the outermost ejecta and the nearby CSM/ISM. The X-ray 
radiation from these shocked layers and the cosmic-rays accel- 
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Fig. 7. Total 511 keV lightcurve for the two remnants of core- 
collapse SNe, Cas A and SN1987A, combining the contribution 
of ^^Ni and '*'*Ti positrons. Also shown are the upper-limits on 
the 5 1 1 keV flux from INTEGRAL/SPl. 



erated at the shock fronts may ionise the freely-expanding ejecta 
to a certain degree. If most of the remnant is ionised, the de- 
cay positrons would experience a more eflicient slowing-down 
(by Coulomb losses), but their annihilation timescale may also 
be longer (if positronium formation in-flight cannot efficiently 
take place). The shocks and shocked layers are also places of 
magnetic field amplification and turbulence where the positrons 
escaping from the freely-expanding ejecta may be trapped. How 
these effects would modify the predicted signals from positrons 
annihilating in the SN/SNR ejecta is currently unexplored. For 
positrons annihilating in the CSM/ISM, the main uncertainty 
is the ionisation state of the medium. In case of strong ionisa- 
tion out to large distances, the annihilation timescale in a typical 
1 cm"^ medium exceeds the ages of the young SNRs consid- 
ered in this work. In addition, as mentioned above, a fraction of 
the escaping positrons can be trapped near the shock fronts and 
may be reaccelerated there (especially since they are akeady in 
the supra-thermal domain). Such effects would allow higher es- 
cape fractions than the above values, without conflicting with the 
INTEGRAL/SPl observations. 



7. Conclusion 

We simulated the 511 keV lightcurves resulting from the anni- 
hilation of the decay positrons of ^^Ni and '^"^Ti in SNe/SNRs 
and their surrounding CSM/ISM. We used a simple model of 
uniform ejecta in homologous expansion and a uniform density 
medium. Two extreme scenarios of positron transport were con- 
sidered: either the positron is free to escape the ejecta along ra- 
dial trajectories, or it is completely trapped in the ejecta. 

We computed specific 511 keV lightcurves for Cas A, Tycho, 
Kepler, SN1006, Gl. 9+0.3 and SN1987A, and compared these 
to the upper-limits derived from INTEGRAL/SPl observations. 
The predicted 511 keV signals from positrons annihilating in the 
SN/SNR ejecta are in all cases below the sensitivity of the SPI 
instrument by several orders of magnitude and therefore do not 
allow constraining the transport of positrons. From the predicted 
511 keV signals for positrons escaping the ejecta and annihi- 
lating in the CSM/ISM, we obtained upper-limits on the ^''Ni 
positron escape fraction of ~13% for Cas A, ~12% for Tycho, 
-30% for Kepler and -33% for SN1006. The 12% upper Hmit 
on the ^''Ni positron escape fraction in a SNIa is consistent with 
the recently estimated value of 5 ±2% required for a nucleosyn- 
thesis origin of the positrons that give rise to the diffuse Galactic 
emission at 5 1 1 keV. 
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